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Abstract 

Two methods of the luminosity determination for the experiment performed by WAS A col¬ 
laboration to search for ^He-?7 bound state are presented. During the measurement the tech¬ 
nique of continous change of the beam momentum in one accelerator cycle (called ramped 
beam) was applied. This imposes the requirement to determine not only the total integrated 
luminosity, but also its variation as a function of the beam momentum. 


1 Introduction 


The existence of 77 -mesic nuclei in which the 77 meson is bound within a nucleus via the strong 
interaction was postulated in 1986 by Haider and Liu [^. Since then 77 - and 77 '-mesic bound states 
have been searched for in many laboratories [2}jT7]. Recent theoretical investigations e.g. 


24 support the search for 77 and 77 '-mesic bound states, however, so far no firm experimental 
confirmation of the existence of mesic nuclei has been found. The discovery of this new kind 
of an exotic nuclear matter would be very important for better understanding of the 77 and 77 ' 
meson properties and their interaction with nucleons inside nuclear matter [^. Furthermore 
it would provide information about the A* (1535) resonance 26 , as well as about the flavour 


singlet component of the quark-gluon wave function of the 77 and 77 ' mesons . 

In November 2010 the search for the ^He -77 bound state was performed with WASA-at-COSY 
facility 28 by measuring the excitation functions for dd and dd ^Hep 7 r“ reactions 


near 


the 77 production threshold [2 [2^^. The measurement was carried out with a deuteron 


beam momentum ramping from 2.127 GeV/c to 2.422 GeV/c, corresponding to the range of the 
excess energy Qg(- 70,30) MeV. During an acceleration process the luminosity could vary due to 
beam losses caused by the interaction with the target and with the rest gas in the accelerator 
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beam line, as well as due to the changes in the beam-target overlap correlated with momentum 
variation and adiabatic shrinking of the beamsize. Therefore, it is necessary to determine not 
only the total integrated luminosity but also its dependence on the excess energy. 

The total integrated luminosity is determined based on the dd ^Hen and quasi free pp pp 
reactions for which the cross sections were already experimentally established. Because of the 
acceptance variation for the beam momentum range for which ^He ions are stopped between two 
Forward Detector layers, the excess energy dependence of the luminosity is determined based on 
quasi-free pp pp reaction for which the WASA acceptance is a smooth function of the beam 
momentum. 

In this contribution we present the procedure of the calculation of the integrated luminosity and 
the determination of the luminosity dependence of the excess energy. 


Determination of luminosity based on the dd —^^Hen 
reaction 


The absolute value of the integrated luminosity was determined using the experimental data on 
the dd ^Hen cross-sections measured by SATURNE collaboration for four beam momenta 
in range between 1.65 and 2.49 GeV/c [^. The cross section (Jdd^^Hen depende nce on the 
transferred momentum squared t = (Ps^g — P^g^ 


may be parametrized as follows 34 35 


d(7{t tmax) 

dt 


= E 


.p^id trnax^ 


die 


(1) 


i=l 


where parameters and bi are described as a function of the total energy y/Sdd'- 


Pdri{y/s^) 


Pi 

\l ^dd Qi 


( 2 ) 


where the values of and were determined by the fit of the above formula to the cross 
sections measured at SATURNE [M| . 

Based on the above parametrization we may determine angular dependence of the cross section 
using a following relation: 


da da dt 

d{cosO'^) dt d{cosO'^) 

where the Jacobian term = 2 • Ip^g^^l • IPs^J is calculated based on the transferred 

momentum squared in the CM system: 


t = {FiHe - ^beamf = “ 2 ■ + 2 ■ \pCeam\ ’ ' COsO*, (4) 

where is the ^He emission angle in the CM frame. 

The available experimental data closest to the range of beam momentum used in the WASA-at- 
COSY experiment for the angular range relevant for our analysis are shown in Eig. Superim¬ 
posed lines present results of the above described parametrisations for beam momenta the same 
as experimental points (red and black) and for two exemplary momenta corresponding to Q=0 
and Q=-40 MeV. 
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Figure 1: (left) Differential cross section as a function of cos^* for SATURNE experimen¬ 
tal data (squares/red and dots/black points for fixed beam momentum p^eam=l-992 GeV/c 
and p^eam=2.492 GeV/c, respectively) and obtained from parametrization (top solid/red, 
dashed/green, dash-dotted/blue, and bottom solid/black lines for Pbeam equal to 1.992 GeV/c, 
2.217 GeV/c, 2.335 GeV/c and 2.492 GeV/c, respectively), (right) The ratio — / d(cose*) 

for p^eam=l-992 GcV/c (squares/red) and p&eam=2.492 GeV/c (dots/black) fitted with second 
degree polynomial functions (dashed/red and solid/black lines, respectively). The marked errors 
result from the statistical experimental uncertainties. 


In the angular region of interest the experimental points he below the curves. Therefore, the 
correction was applied for the ^He angular range from about 4° to 10° which corresponds to the 
cos^* G(0.88,0.98) for the considered reaction. The ratio between experimental (SATURNE) and 
parametrized cross section — / d(cos9*) fitted with second degree polynomial function 

for both experimental beam momentum values: 1.992 GeV/c and 2.492 GeV/c. Obtained result 
is presented in the right panel of Eig. The cross section correction A is calculated for fixed 
cos^* using the fitted functions and interpolated for the proper beam momentum value from 
range Pbeam ^ (2.127, 2.422)GeV/c. 

The measurement of the dd Hen reaction was based on the registration of the outgoing helium 
in the Eorward Detector. Low-energetic ^He ions were stopped in the 3rd layer of the Eorward 
Range Hodoscope, while high-energetic ions were stopped in the 4th layer. The helium identifica¬ 
tion was based on the AE-AE method. The outgoing neutrons were identified using the missing 
mass technique. In order to reduce background originating from quasi-free dp{n) HenTr^, the 
cut in missing mass rrix vs. missing energy Ex spectrum was applied as it is presented in upper 
panel of Eig. Additionally, for high beam momentum region background was subtracted via 
fitting the signal and background function to the missing mass spectrum for different intervals 
of cosO* and beam momentum, what is presented in the lower panel of Eig. 

In order to calculate the total integrated luminosity, the number of events, efficiency, as well as 
cross section was determined for 5 intervals of cosO* in the range from 0.88 to 0.98 and 5 intervals 
of excess energy Q in the range from -70 MeV to 30 MeV corresponding to the angular range 
of the reaction and the beam momentum ramping, respectively. The integrated luminosity was 
then calculated for each (i,jf)-th interval in following way: 


T int _ 




Ni,, 


d(j i j j 

d(cos9 *) 


A(cos^*) 


(5) 
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Figure 2: (upper panel) The missing mass rrix vs. missing energy spectrum for simulations 
(left) and DATA (right). Applied cut is marked with red line, (lower panel) The missing mass 
rrix spectrum for i. e. cos6>* G(0.96,0.98) and Qg( 0,5) MeV. The red line shows fit to the signal 
and background while green line shows fit of the Gauss function to the background. Signal peak 
is marked as a blue line. The main background on the right side of the black line corresponds to 
the quasi-free dp{n) HeriTr^ reaction. 


where A{cos0*) is the width of the cos^* interval. The overall efficiency including reconstruction 
efficiency and geometrical acceptance of the detector was determined based on the Monte Carlo 
simulations and is varying between 50% and 70%. 

The preliminary luminosity dependence of co5^* for whole excess energy range is presented in 
Fig-i The total integrated luminosity was calculated as a weighted average of the luminosities 
determined for individual cosO"^ intervals: 


T tot _ 

^dd^^Hen ~ 


l^i=l {ALj) 

E 5 1 

i=l (AL^)2 


AT 


tot _ 

dd^^Hen ~ 



- 1/2 


(6) 


The average integrated luminosity with its statistical uncertainty equals L^f^^3^^^={1102E2)nb ^ 
It is marked in Fig. with dashed red line. 
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Figure 3: Integrated luminosity as a function of cosO"^. The statistical uncertainties are marked as 
a vertical bars. The preliminary established weighted average of integrated luminosity is marked 
as a dashed red line and is equal to 1102±2n6“^ where only a statistical error is given. The 
analysis was carried out with condition that the number of ’’neutral clusters” in reconstructed 
in the Central Detector is less than 2. 


3 Luminosity dependence on the excess energy 


In order to determine the luminosity dependence on the beam momentum we used the quasi¬ 
elastic proton-proton scattering in the deuteron-deuteron collisions: dd ppUspUsp. In this 
reaction protons from the deuteron beam are scattered on the protons in the deuteron target. 
We assume that the neutrons are acting only as spectators which means that they do not take 
part in reactions but move with the Fermi momentum of their parent deuterons. 

In the case of quasi-free proton-proton scattering the formula for the calculation of the integrated 


luminosity can be written in the following form 36 : 


L = 


NoN,,^ 




dcr ( 

An(eiab,(t^lab) 


0*, Pf^o , o : ^Fio)f {PFio 5 ^Fio 5 ^Fio )dpFrodcOs0F^od(l)Fx- 


, dip'^dcosO* 

(7) 


The formula is determined based on the fact, that the number of quasi-free scattered protons 
into the solid angle ACl{0iab^ (piab) is proportional to the integrated luminosity L, as well as the 
inner product of the differential cross section for scattering into the solid angle around 6>* and 0* 
angles expressed in proton-proton CM system: 5 5 ^^^ 1 , 2 ) probability 

density of the Fermi momentum distributions: f{pFi 2 ? 25 2 ) iiiside the deuteron beam and 

deuteron target, respectively. The detailed description of the luminosity calculation for quasi-free 
reaction one can find in Ref. [^ . 

Due to the complex detection geometry a solid angle corresponding to particular part of the 
detector cannot be in general expressed in a closed analytical form. Therefore, the integral 
in above equation was computed with the Monte-Carlo simulation programme, containing the 
geometry of WASA detection system and taking into account detection and reconstruction effi¬ 
ciencies. The Monte-Carlo simulations were carried out for the deuteron beam momentum range 
Pbeam G (2.127,2.422)GcV/c corresponding to the experimental ramping. The program first choose 
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randomly the momentum of the nucleon inside the deuteron beam and deuteron target, respec¬ 
tively, according to the Fermi distribution 
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Then, the total proton-proton invariant mass 
yopp and the vector of the center-of-mass velocity are determined. Next, the effective proton 
beam momentum was calculated in the frame where one of the proton is at rest and mo¬ 

mentum of protons is generated isotropically in the proton-proton center-of-mass frame. Further 
on, the momenta of outgoing particles are transformed to the laboratory frame and are used 
as an input in the simulation of the detection system response with the GEANT computing 
package. For each of Nq simulated event we assign a weight corresponding to the differential 
cross section, which is uniquely determined by the scattering angle and the total proton-proton 
collision energy s/s^. 

The factor Nq/2t: in Eq. [^is a normalization constant. It results from the fact that the integral 
is not dimensionless and its units correspond to the units of the cross sections used for the calcu¬ 
lations. Therefore, it must be normalized in such a way that the integral over the full solid angle 
equals to the total cross section for the elastic scattering averaged over the distribution of the 
total proton-proton invariant mass resulting from the Eermi distribution of the target and 


beam nucleons. In the absence of the Eermi motion it should be simply equal to a total elastic 
cross section for a given proton beam momentum. A factor 27r comes from the fact that protons 
taking part in the scattering are indistinguishable. 


The differential cross section for quasi free dd pp^spUsp reaction is a function of the scattering 
angle 0* and the total energy in the proton-proton centre-of-mass system which is dependent 
on effective proton beam momentum pH^I^ seen from the proton in the proton-proton system. 
In order to calculate it, we have used the cross section values for proton-proton elastic scattering 
pp pp computed based on the SAID program 38 because the EDDA collaboration data 
base 39 is insufficient. The distribution of the effective beam momentum as well as a comparison 


of the SAID calculations and the existing differential cross section from the EDDA measurements 
are shown in Eig. As we can see, the differential cross sections calculated using the SAID 
programme are in agreement with distributions measured by the EDDA collaboration. 




Eigure 4: (left) Differential cross sections for proton-proton elastic scattering as a function of the 
beam momentum for a three values of the scattering angle 0* in the CM frame. Black points 
show EDDA collaboration data [^, while lines denote SAID calculations [^. Distribution of 
the effective beam momentum for quasi-free pp pp reaction calculated for the deuteron beam 
momentum range Pbeam G(2.127,2.422)GeV/c is also presented in the figure, (right) Bilinear 
interpolation of the differential cross section figure is adapted from 


36 
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The differential cross section for appropriate and 6>* was calculated using bilinear interpo¬ 

lation in the momentum-scattering angle plane according to the formula: 


= (1 - i)(l - 


0*’^) +t{l - + 


dVt 


tu^^iPbeam^O*’^) + {1 - t)u^{pl^ 


(8) 


,^’ 2 ) 


where t and u variables are defined in right panel of Fig. 

The number of experimental events N^xp was determined based on conditions and cuts described 
in details in reference 40 . In the analysis, at the beginning, we carried out primary events 


selection applying condition of exactly one charged particle in the Forward Detector (FD) and 
one particle in the Central Detector (CD). 

In Ref. we can find detailed studies of the possible background reaction contributions. In 
case of this analysis the dominating background processes are dd pprispUsp dTr+ngpUsp, 
dd diyPtfisp and dd PPsp'^'^sp^ where the subscripts sp, b and t denote the spectators, particles 
from the beam and from the target, respectively. In order to reject the events corresponding 
to the charged pious registered in the Central Detector, the cut on the energy deposited in the 
Electromagnetic Calorimeter (Cal) vs. energy deposited in Plastic Scintillator Barrel (PSB) 
spectrum was applied and is presented in the Fig. 



Edep(PSB) [GeV] 


Figure 5: Experimental spectrum of the energy loss in the Plastic Scintillator Barrel shown as a 
function of the energy deposited in the Electromagnetic Calorimeter. The applied cut is shown 
as a black line. Pious in data spectrum are concentrated for Edep(PSB) around O.OOSGeV. 

It is not possible to separate quasi-elastic p-p scattering from the quasi-elastic d-p scattering, 
however it was investigated that for the forward scattering angles of about the d-p 

cross sections are about 20 times smaller than p-p cross sections and we take this uncertainty of 
abut 5% as a systematic error to the final result. The applied cut in polar angle OpD is shown 
in Eig. In order to subtract the background coming from dd Pbdt'^sp reaction, the range 
dcD G (40,100)° was taken into account in further analysis. 

Additionally, the background was subtracted in A0 = (t>FD—(t>CD spectrum. In order to symetrize 
the background instead of | A0| we define (27r-hA(/))mod27r. Afterwards, the background was fitted 
with 1st order polynomial for each of excess energy Q intervals. The exemplary (27r + A0)mod27r 
spectrum is presented in Eig. 
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Figure 6: Correlations between the polar angles OpD and Ocd for the WMC Simulations of 
dd pprispUsp reaction (left panel) and obtained in experiment (right panel). Applied cut is 
marked with red line. The indicated area correspond to the: a) dd ppUspUsp, b) dd d^ptfisp 
and dd pprigpUsp^ c) dd pPsp'^ingp^ d) dd PbdtTisp- 



( 27 r+A(|))nnod 27 r [deg] 


Figure 7: Distributions of (27r+A(/))mod27r, where A0 = (j)cD — (t>FD is the difference of azimuthal 
angles in Central Detector and Forward Detector. The example spectrum for one of the Q intervals 
(black line) with marked fit function (red line) and signal peak after background subtraction (blue 
line) is presented. 


After all cuts and conditions described above, the number of experimental data was determined 
and the luminosity was calculated according to formula 0 for each excess energy interval taking 
into account also prescaling factor of the applied experimental trigger equal to 4000 as well as 
shadowing effect equal to 9%. The latter results from the fact that proton is shadowed by the 
neutron inside the deuteron which reduces the probability of the quasi-elastic scattering. Un¬ 
fortunately, there are no experimental results about the shadowing in dd pprigpUsp collisions. 
However, we can try to estimate it based on the probability that a neutron shadows the proton 
in one deuteron which equals 0.045 and assume that shadowing appears independently in 
deuteron beam and deuteron target. The rough estimation of the probability that the shadowing 
will not take place in dd reaction (1 - 0.045)^ gives about 0.91. 

The preliminary result is presented in Fig. The statistical uncertainty of each point is about 
1%. The luminosity variation (increase in the excess energy range from about -70MeV to - 
40MeV, and then decrease) is caused by the change of the beam-target overlapping during the 
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acceleration cycle and also by adiabatic beam size shrinking 42 . The obtained total integrated 
luminosity within its statistical uncertainty is equal to = . For further 

analysis the luminosity was fitted by 3 degree polynomial + 6(5^ TcQ + d. The fitted function 

is marked with the red line in Fig. 
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Figure 8: Integrated luminosity calculated for experimental data for quasi-free dd pp^spUsp 
reaction (blue points) with fitted 3 degree polynomial function (red line). 


4 Systematics 


In case of the dd ^Hen reaction one source of the systematic error originates from the variation 
of the cuts used for separation of high-energetic helium in Forward Detector and is equal to about 
2%. Additionally we have also taken into account an uncertainty due to the method used for the 
background subtraction amounting to 1.6%. Another source of the luminosity calculation error 
is connected to normalization to SATURNE experiment and originates from three independent 
sources: i) statistical error of the SATURNE data (6.5%), ii) normalization uncertainty of the 
SATURNE data for the dd ^Hen cross sections (7%) and iii) assumption of linear interpolation 
between SATURNE points used for the estimation of the correction for the parametrized cross 
section presented in Eig. (<1.8%). 

The systematical errors for dd pprispUsp analysis resulting from the change of the cuts used for 
the separation of the quasi-free pp scattering from the background (Eig. [^and Eig. is equal to 
to about 4.1%. The another contribution to the systematical error comes from the assumption of 
the potential model of the nucleon bound inside the deuteron and is equal to about 0.8%. This 
uncertainty was established as the difference between results determined using the Paris 
and the CDBonn 43 potentials. The next source of the systematic error may be attached to 


the assumption of the bilinear approximation of the cross section shown in Eig. (right). This 
systematical uncertainty was estimated using assumption in which instead of the interpolation 
we took the cross section value from the closest data point in the effective proton beam mo¬ 
mentum - scattering angle plane. The performed calculations give the difference of about 1.8%. 
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Additionally we have also taken into account an uncertainty related to the background subtrac¬ 
tion in (27r + A0)niod27r spectra which does not exceed 0.6%. The systematical uncertainty 
includes also contribution connected to the shadowing effect. Untill now, we have no theoretical 
estimation of the possible error of this effect calculation, therefore conservatively we take as an 
systematic uncetainty half of this effect: 4.5%. In the systematical error calculation we also take 
into account the uncertainty 5% resulting from the background of the quasi-elastic d-p scattering 
(Sec. [^. The normalization error includes two contributions: normalization uncertainty of the 
EDDA data (4%) and the systematical errors for pp elastic scattering cross-sections (2.7%) [^ . 

The cross section was approximated by the calculation using the SAID procedure. Because, the 
SAID cross section very well describes EDDA data, we assume the systematical errors of the 
differential cross section based on EDDA calculations. 

The total integrated luminosity calculated based on dd ^Hen and the quasi-free dd pprispUsp 
reactions with statistical, systematical and normalization error are equal to = (1102 ± 

^stat ‘^^syst lO^T^orm)'^^ and = (1329 dz ‘Igtat ^ ^^^syst ^ i respec¬ 

tively. The systematical and normalization errors were calculated by adding in quadrature the 
appropriate contributions described above. 

5 Summary 

We carried out the luminosity determination for the experiment performed with WASA-at-COSY 
to search for the ^He-r^ bound states in deuteron-deuteron fusion. The luminosity was calculated 
based on two reactions: dd ^Hen and the quasi-free dd pprispUsp. We calculated the total 
average integrated luminosity and compared it for both channels. The obtained results are 
consistent, however within large normalization errors. 
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